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Abstract. Using the Very Large Telescope in Multi Object Spectroscopy mode, we have observed a sample of 1 13 field spiral 
galaxies in the FORS Deep Field (FDF) with redshifts in the range 0. 1 < z < 1 .0. The galaxies were selected based on apparent 
brightness (R < 23"') and encompass all late spectrophotometric types from Sa to Sdm/Im. Spatially resolved rotation curves 
have been extracted for 77 galaxies and fitted with synthetic velocity fields taking into account all observational effects from 
inclination and slit misalignment to seeing and slit width. We also compared different shapes for the intrinsic rotation curve. 
To obtain robust values of V milx , our analysis is focused on galaxies with rotation curves that extend well into the region of 
constant rotation velocity at large radii. If the slope of the local Tully-Fisher relation (TFR) is held fixed, we find evidence 
for a mass-dependent luminosity evolution which is as large as up to AM B x -2"' for the lowest-mass galaxies, but is small 
or even negligible for the highest-mass systems in our sample. In effect, the TFR slope is shallower at z ~ 0.5 in comparison 
to the local sample. We argue for a mass-dependent evolution of the mass-to-light ratio. An additional population of blue, 
low-mass spirals does not seem a very appealing explanation. The flatter tilt we find for the distant TFR is in contradiction to 
the predictions of recent semi-analytic simulations. 
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1. Introduction 

Ever since the relation between the luminosity L and the max- 
imum rotation velocity V max of spiral galaxies was first ob- 
served (Tully & Fisher IT977> . the physical origin of its slope 
and scatter, as well as the possible evolution thereof over dif- 
ferent cosmic epochs have been subject to debate both in theo- 
retical and observational studies. Within the last few years, the 
Tully-Fisher relation (TFR) has been put into the framework 
of a Fundamental Plane (FP) for spiral galaxies that introduces 
the disk scale length as a third parameter (e.g. Burstein et 
al.Qp97i. Similiar to the FP of dynamically hot galaxies, i.e. 
stellar systems that are stabilized due to random motion (e.g. 
Dressier et al. 1987), the spiral FP has smaller scatter edge-on 
than its two-dimensional projections. This may be understood 
in terms of two dominant parameters for disk galaxies, for ex- 
ample the mass and angular momentum (Koda, Sofue & Wada 
20001. 
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Numerical simulations within hierarchical Cold Dark 
Matter (CDM)-dominated cosmologies have been successfully 
used to reproduce the observed slope of the local TFR, whereas 
the numerical zero points were offset due to dark halos with too 
high concentrations (e.g. Navarro & Steinmetz 2000). The TFR 
slope is predicted to remain constant with cosmic look-back 
time in such N-body simulations; nevertheless the modelling 
of realistic stellar populations at sufficient resolution remains a 
challenge. 

Other theoretical approaches focus more on the chemo- 
spectrophotometric aspects of disk galaxy evolution. For ex- 
ample, Boissier & Prantzos (2001 ) used the "hybrid" approach 
(Jimenez et al. 1998 1 that relates the disk surface density to 
the properties of the associated DM Halo, and calibrated it to 
reproduce the observed colors of local spirals. Compared to 
these, the authors predict higher luminosities for large disks 
and lower luminosites for small disks at redshifts z > 0.4. A 
similiar evolution is found by Ferreras & Silk ( 2001 1. By mod- 
elling the mass-dependent chemical enrichment history of disk 
galaxies with the local TFR as a constraint, the authors find a 
TFR slope that increases with look-back time (i.e., for a pa- 
rameter! sation L oc V^ ax , a increases with redshift). 

In the last decade, many observational studies of the lo- 
cal TFR have produced very large samples with N \,j ~ 1000 
(e.g. Haynes et al. 1999 1, not only to derive the slope and scat- 
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ter with high accuracy, but also to map the peculiar velocity 
field out to cz ~ 15000km s _1 (e.g. Mathewson & Ford 1996 1. 
Other groups used spirals, partly with cepheid-calibrated dis- 
tances, to measure the Hubble constant. For example, Sakai et 
al. (|5000l derived a value of H = (71 + 4)kms~ I with this 
method. 

At higher redshifts, robust measurements of rotation ve- 
locities and luminosities become increasingly difficult. This is 
partly because of the low apparent magnitudes of the galax- 
ies, but also due to the limited intrinsic spatial resolution (see 
Sect.0]for a more detailed description of this effect). A number 
of samples with 10-20 objects in the regime 0.25 < (z) < 0.5 
have been observed in recent years to estimate a possible evo- 
lution in luminosity by comparison to the local TFR. The re- 
sults of these studies were quite discrepant: e.g. Vogt et al. 
( 1996| 119971 find only a modest increase in luminosity of 
AM B « -0.5 m , whereas Simard & Pritchet ( 1998 1 and Rix et al. 
dl997> derive a much stronger brightening with AMg ~ -2.0™. 
A study of 19 field spirals by Milvang-Jensen et al. (2003 1 re- 
veals a value of AMb ~ -0.5'" and shows evidence for an in- 
crease with redshift. Another sample of 19 spirals by Barden et 
al. (|2003 1 which covers the high redshifts 0.6 < z < 1 .5 yields 
a value of AM B « -l.l m . 

It seems likely that some of these results are affected by the 
selection criteria. For example, Rix et al. selected blue colors 
with (B-R) ohs < 1 .2'", Simard & Pritchet strong [O II] emission 
with equivalent widths >20A, while Vogt et al. partly chose 
large disks with r^> 3 kpc. The two former criteria prefer late- 
type spirals, whereas the latter criterion leads to the overrepre- 
sentation of large, early-type spirals. Additionally, due to the 
small samples, all these studies had to assume that the local 
TFR slope holds valid at intermediate redshift. We will further 
discuss this issue in Sect.[S] 

Based on a larger data set from the DEEP Groth Strip 
Survey (Koo 2001 1 with N ~ 100 spirals in the range 
0.2 < z < 1.3, Vogt (|2001 1 finds a constant TFR slope and a 
negligible rest-frame fi-band brightening of less than 0.2 mag. 
In a more recent publication from this group which investigates 
the luminosity-metallicity relation, an evolution both in slope 
and zero point is observed (Kobulnicky et al. 2003 1, in the sense 
that the luminosity offsets are largest at the low-luminosity end 
of the sample and smallest at the high-luminosity end. The au- 
thors argue that low-luminosity galaxies could have either un- 
dergone a decrease in luminosity or an increase in the metallic - 
ity in the last ~ 8 Gyrs. 

Preliminary results from our TF project indicating a mass- 
dependent luminosity evolution of distant field spirals have 
been presented in a letter (Ziegler et al. 2002 1. In this paper, we 
will describe the derivation of the maximum rotation velocities, 
the galaxies' structural parameters and the luminosities in more 
detail and present the data table of the enlarged, full sample. 
Complementary to our approach in Ziegler et al., the complete 
analysis will be restricted here to galaxies with rotation curves 
that extend well into the region of constant rotation velocity 
at large radii, i.e. spirals that yield robust values of V max . In 
addition, different shapes for the intrinsic rotation curves will 
be compared. Finally, we will discuss potential environmental 
effects on the sample and present the galaxies' virial masses. 



The paper is organized as follows. In Sect. |3 we describe 
our selection procedure and the observations. Data reduction 
will be outlined in Sect. [5] The extraction of rotation curves 
and V max derivations are described in Sect. 0] followed by the 
details of the transformations from apparent to absolute magni- 
tudes in Sect.|5]and the presentation of the data table in Sect. [6] 
We will then construct the distant TFR in Sect. and discuss 
the results in Sect. [8] A summary is given in Sect. [9] 

Throughout this article, we will assume the concor- 
dance cosmology with Q. m = 0.3, £2a = 0.7 and Hq 
70 kin s 1 Mpc -1 . 

2. Sample selection and observations 

Our sample consists of galaxies in the FORS Deep Field (FDF), 
a sky region near the south galactic pole with deep UBgRIJK 
photometry and visible completeness limits similiar to the 
Hubble Deep Fields. For a description of the field selection cri- 
teria and the total N « 8750 object catalogue of the FDF, we 
refer to Heidt et al. (2001}. 

The basis for spectroscopy target pre-selection was the 
FDF photometric redshifts catalogue (Bender et al. 2001). To 
keep the selection function as simple as possible, the only spec- 
trophotometric costraints were that the galaxies should have a 
spectral energy distribution (SED) later than E/S0 and an appar- 
ent total brightness of R < 23™. This limit was chosen to gain 
S/N « 5 in emission lines at intermediate resolution within 2-3 
hours integration time with the Very Large Telescope (VLT). 
To ensure the visibility of either the [O II] 3727 doublet, HB 
or [O III] 5007 within the wavelength range of the 600R grism 
of FORS, the upper limit for the photometric redshift was 
Zphot ^ 1-2. Since the basic aim was to derive spatially resolved 
rotation curves of the galaxies, (most of the) objects with incli- 
nations i < 40° were rejected. 

In Multiobject Spectroscopy (MOS) mode, FORS offers 19 
individually moveable slits. The setups were prepared with the 
FIMS (FORS Instrument Mask Simulator) package. According 
to the position angles, we subdivided the target galaxy sample 
into bins of 30° to minimize geometric distortions in the final 
rotation curves. Thus, a set of 6 MOS masks was necessary to 
cover all orientations. In the case of some objects, either the 
inclination limit of 40°, the limit of 15° deviation between po- 
sition angle and slit orientation, or the magnitude limit had to 
be exceeded to fill all slitlets. 

The first observations were carried out with FORS2 
mounted on VLT Unit Telescope 2 (UT2) in September and 
October 2000. Each mask was exposed for 3 x 3000 s. The slit 
widths were set to one arcsecond and grism 600R with order 
separation filter GG435 was used, yielding a spectral resolu- 
tion of R S3 1200. FORS was operated in its standard reso- 
lution configuration of 0.2"/pixel. The seeing conditions cov- 
ered the range 0.43" < FWHM < 0.81" according to 
the Differential Image Motion Monitor (DIMM). In October 
2001, an additional set of three MOS masks was observed with 
FORS 1 mounted on VLT UT3 with the same instrument con- 
figuration as in 2000 with seeing of 0.74" < FWHM < 0.89". 
In total, spectra of 129 spirals out of 156 candidates were taken 
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(the latter number includes targets which slightly exceed the 
brightness limit). 

3. Data reduction 

The reduction procedures were implemented in the ESO- 
MIDAS (Munich Image Data Analysis System) environment. 
For all three observing runs, the bias sets taken at the beginning 
and end of the night showed a very stable two-dimensional 
structure, so all bias frames from the nights of the 2000 and 
2001 observations were used to generate two master biases. 
The individual bias frames have been normalized to the same 
median count rate and were median-averaged, followed by a 
slight Gaussian filter smoothing (<x r v = 3 pixel). The overscan 
region of the science frames was used to determine the offset 
constant for final master bias subtraction. Since the dark current 
was very uniform across the CCDs, an extra dark subtraction 
has been neglected. 

Spectroscopy flatfield (FF) through-mask exposures with 
FORS were taken by default with two sets of lamps switched 
on alternately. Since either the upper or lower half of the frames 
were not usable due to contaminating light from the gaps be- 
tween the slits, the appropriate regions from both sets of FF 
frames were extracted and re-combined after multiplicative 
normalization and median-averaging. The slit regions of cal- 
ibration and science frames were then individually extracted. 
The FF slit exposures were approximated by a polynomial fit of 
sixth degree in dispersion direction (X axis) to account for the 
CCD response curve, which was used for normalization prior 
to the correction of the pixel-to-pixel variations. 

In the next step, the geometric distortions caused by the fo- 
cal reducer were corrected. For the slits at the bottom and top 
of the CCD, where the distortions were maximum, the curva- 
ture of spectral features corresponded to a displacement of up 
to 5-6 pixel both in X and Y direction (see Fig.^for an illustra- 
tion). A polynomial fit of second degree was fitted to the galaxy 
spectra in each slit to derive the curvature along the spatial axis. 
Based on this fit, the science and wavelength calibration spec- 
tra were rectified with an accuracy of 0. 1 pixel (corresponding 
to 0.02"). Tests revealed that the flux conservation of this pro- 
cedure was accurate to within a few percent. The distortions 
along the X axis could be corrected during standard wavelength 
calibration. In the calibration exposures, the HgCd lamp was 
switched on additionally to the He, Ar and Ne lamps in order 
to gain a sufficient number of emission lines below ~ 5800 A. 
For the two-dimensional dispersion relation, polynomial fits of 
third and first degree were used in the direction of the X and Y 
axes, respectively. The typical r.m.s. of the relation was 0.03- 
0.04 A at a stepsize of 1 .08 A per pixel. 

The night sky emission was fitted column by column with 
first order fits, unless a galaxy spectrum was located at the ex- 
treme edges of the slit. In those cases, zero order fits yielded 
the best results. Night sky subtraction was performed individ- 
ually for each mask exposure. Prior to the final addition of the 
three exposures, the optical center of each galaxy along the Y 
axis has been determined by fitting a Gaussian. If necessary, the 
spectra were shifted by integer values (two pixels at maximum 
in a few cases) to ensure consistent profile centers to within at 




Fig. 1. Upper part of a raw FORS 1 MOS frame with 3000 sec- 
onds exposure time showing the region of slit #2. Wavelength 
is increasing parallel to the X axis from left to right. A rela- 
tively bright (R = 19.94'") galaxy spectrum is overlaying the 
night sky emission lines. Note that the magnification scale is 
larger along the Y axis to demonstrate the distortions by the 
focal reducer in this direction. 



least half a pixel, corresponding to 0.1". A weighted addition 
was used if the average seeing varied by more than ~ 20 percent 
between the exposures. 

4. Rotation velocity derivation 

4.1. Redshift distribution 

Out of the 129 galaxies of which spectra were taken, red- 
shifts could be determined for 113 spirals, including three cases 
of secondary objects which were covered by a MOS slit by 
chance. The 16 galaxies without spectroscopic redshifts popu- 
late the extreme faint end of our apparent brightness distribu- 
tion with a median of (R) = 22.9"'. Since only one of these 
spirals features colors of a very late-type SED, the S/N of the 
remainders may be just too low to yield detectable emission 
lines. According to the photometric redshifts, the [O II] 3727 
doublet is possibly redshifted out of the wavelength range of 
the R600 grism in the case of 8 targets, while for 5 other galax- 
ies in the regime 0.2 < z p hot < 0.5, only [OIII] or H/3 emission 
is potentielly covered, which is weaker than Ha or [Oil] for 
typical spirals. 

Fig. |3 shows the redshift distribution of our sample, re- 
stricted to objects with appropriate rotation curves for the TF 
analysis, see Sect. 14. 41 for the constraints. The median redshift 
is (z) = 0.45. It is likely that the bimodal shape is a combina- 
tion of the two following effects. 

Firstly, as has been outlined in Heidt et al. (2003 1, the 
southwestern corner of the FDF covers the outskirts of a galaxy 
cluster at z = 0.33. Allowing a spread in redshift of Az = 0.01, 
we find that a maximum of 8 galaxies in our sample could be 
members of this cluster. The small redshift "bump" at z ~ 0.3 
can be attributed to these objects. 
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We therefore also performed a comparison to the SED tem- 
plates in the Kennicutt catalogue ( 1992a) with a focus on the 
relative line strengths of the Hf3/[0 III] emission. As a third in- 
dicator, the rest-frame [O II] 3727 equivalent widths were de- 
rived and correlated to type with the values for local spirals 
given by Kennicutt (1992b I as a reference. 

A combination of these three criteria yielded 14 spirals of 
type T < 3 in the redshift range 0.11 < z < 0.89 with 
a median of (z) = 0.42, 43 galaxies with T — 5 covering 
0.09 < z < 0.97 with (z) = 0.45 and 20 objects with T > 8 
in the regime 0.23 < z < 0.97 with (z) = 0.51. 

4.3. Rotation curve extraction 



Fig. 2. Redshift distribution of the 77 FDF spirals with usable 
rotation curves (thick solid line). The shaded region gives the 
distribution after exclusion of 8 possible members of a cluster 
at z = 0.33. The dashed line represents the lower resolution 
spectroscopy sample observed with grism 1501 as part of the 
FDF high-z study (Noll et al. I2004> . Both data sets feature a 
redshift "gap" around z ~ 0.6, see text for details. 



Secondly, the sample contains relatively few galaxies 
around redshift z ~ 0.6. For comparison, we also show the 
distribution of 144 galaxies at z < 1-1 from the FDF high-z 
campaign (Mehlert et al. 120021 Noll et al. l2004t which are not 
contained in our sample. The high-z spectra were taken with 
the same instrument configuration as the TF data, except that 
the low-resolution grism 1501 was used which covers a much 
broader wavelength range of ~ 6000 A in observer's frame. The 
large number of galaxies with z < 1 which did not enter our TF 
survey can be attributed to the much fainter brightness limits of 
the high-z study (/ < 24.5'"), the inclusion of elliptical galaxies 
and the lack of constraints to the inclinations. It is clear from 
Fig.|3that both samples feature the same redshift "gap". 

As a test, we verified our selection criteria on the most re- 
cent version of the FDF photometric redshifts catalogue, find- 
ing that only a handfull of types later than Sa could have been 
missed by the original selection. We therefore conclude that the 
volume probed with the FDF probably contains fewer galaxies 
at z * 0.6 than the neighbouring redshift bins and that the dis- 
tribution of the TF spirals is unlikely to reflect a selection bias 
or an observational effect. 



4.2. Spectrophotometric classification 

To gain a classification of our spectra, three criteria were used. 
The SED model parameter from the photometric redshifts cat- 
alogue, which is related to the star formation e-folding time of 
the fitted templates, can be transformed into the de Vaucouleurs 
scheme (ranging from T — 1 for Sa to T — 10 for Im). This 
is useful especially for spectra with very few identifiable lines 
and/or low S/N, but has the disadvantage that dust reddening 
might induce a classification of too early a type for highly in- 
clined spirals. 



The measurements of the rotation velocity as a function of ra- 
dius were performed in a semi-automatic manner. In the first 
step, 100 columns of the spectrum centered on the considered 
emission line were averaged to get a profile of line plus con- 
tinuum along the spatial axis. This profile was approximated 
with a Gaussian to derive the optical center to within 0.1". 
For a few objects located very close to the slit edge, the cen- 
ter had to be redefined manually. Then, the emission lines were 
fitted row by row. To enhance the S/N, three neighbouring rows 
were averaged prior to the emission line fitting; for very weak 
lines, this "boxcar" was enlarged to five rows (correspond- 
ing to one arcsec). In the case of the [Oil] 3727 doublet, two 
Gaussians with equal FWHM and an observer's frame separa- 
tion of 2.75 (1 + z) [A] were assumed as a line profile approxi- 
mation, while a single Gaussian was used for [OIII] 5007, H/3 
or Ha, with the latter being visible only in four spectra. The 
red- and blueshifts along the spectral axis due to rotation were 
measured relative to the observed wavelength of the line at the 
optical center and converted into velocity shifts after cosmolog- 
ical correction by a factor of (1 + z) _1 . This position-velocity 
information defines an observed rotation curve (RC). 

Each curve was visually inspected prior to the RC mod- 
elling. Seven objects had too low a S/N to derive spatially re- 
solved rotation velocities. We also rejected RCs with strong 
asymmetries or other signatures of substantial kinematic dis- 
turbances, "solid-body" rotators and objects that did not show 
any rotation at all within the measurement errors. In total, the 
RCs of 77 spirals were appropriate for the V max derivation. We 
present a range of examples in Fig.|5] 



4.4. Rotation curve modelling 

For distant, apparently small galaxies, the effect of the slit 
width on the observed rotation velocity as a function of radius, 
V° b t s (r), must be considered. At redshift z = 0.5, a scale length 
of 3 kpc — typical for an L* spiral — corresponds to approx. 
0.5 arcsec only, which is half the slit width used in our MOS 
observations. Any value of V^ir) is therefore an integration 
perpendicular to the spatial axis (slit direction), a phenomenon 
which is the optical equivalent to "beam smearing" in radio ob- 
servations. If not taken into account, this effect could lead to an 
underestimation of the intrinsic rotation velocities. 
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Fig. 3. Examples of rotation curves from our data. Observed rotation velocity is plotted against the projected distance from 
galaxy center. Error bars denote the errors of the Gaussian fits to the emission lines. Object numbers and emission lines used are 
given for each spiral (either [O n] 3727, [O m] 5007, Ha or H/3). The observed rotation curves have been modelled with synthetic 
velocity fields (solid lines) which take into account geometric distortions due to inclination and slit misalignment as well as 
seeing and optical beam smearing, see Sect. l4.4l for details. The small spatial extent of the curve of galaxy #5515 to the left hand 
side is due to contaminating light from a foreground elliptical. 



We overcame this problem by generating synthetic RCs. 
Our approach is similiar to the procedure described by Simard 
& Pritchet ( 1999 1 but differs from their fitting method in some 
ways. E.g., the scale length of the emitting gas is a fixed pa- 



rameter in our algorithm (see below), and the deviation angle 
between apparent disk major axis and slit direction is taken into 
account. Moreover, in our approach we do not fit simulated 2D 
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spectra to observed 2D spectra but simulated RCs to observed 
RCs. 

For the simulation, one has to assume an intrinsic rota- 
tional law V™[(r). As a first of three variants, we used a simple 
shape with a linear rise of Vjjj! at small radii, turning over into a 
region of constant rotation velocity where the Dark Matter Halo 
dominates the mass distribution. This can be achieved with the 
parameterisation 



v« 



(r a + r ') I/fl 



(1) 



(e.g. Courteau 1997 1 with a factor a that tunes the sharpness 
of the "turnover" at radius r = ro, and a constant rotation 

^r"f( r ) = ^ max f° r r r °- We usec ' a ran 8 e °f values for a 
on a set of ~ 20 RCs and found that a — 5 best reproduced the 
observed shape of the turnover region. To minimize the number 
of free parameters, we kept a fixed to that value for all objects 
in the further analysis. However, due to the heavy blurring of 
the curves, a is the least critical parameter in Eq.^ 

The turnover radius ro was assumed to be equal to the scale 
length of the emitting gas, which is larger than the scale length 
derived from continuum emission. This is discussed in Ryder 



& Dopita ( 19941 and Dopita & Ryder ( 1994i using observa- 
tional and theoretical approaches, respectively. Since the struc- 
tural parameters like the continuum scale lengths were mea- 
sured in an /-band VLT image (see Sect. 15. U . we derived the 
corresponding gas scale length via 



r = (2 -z/2)r d . 



(2) 



This equation yields ro ~ 2 rj for the least distant FDF spirals 
(for which / corresponds to rest-frame R) and ro = 1.5 r& at 
Z = 1, where / corresponds to rest-frame B. In other words, 
Eq. |2 is used to gain a correlation between rest-frame scale 
length and gas scale length that is in compliance with the results 
of Ryder & Dopita. We also tested this equation directly on a 
few high S/N emission line profiles in our spectra. It should be 
noted however, that the V max derivation is much less sensitive 
to the scale length than to the inclination and the misalignment 
angle. This is mainly due to the smoothing by the instrumental 
PSF 

We alternatively used two other templates of intrinsic RC 
shapes within the Universal Rotation Curve (URC) framework 
of Persic & Salucci (119911 cited as the URC91 hereafter) and 
Persic et al. (1996 URC96 hereafter). The authors introduced 
a dependence of the RC morphology on luminosity: Only L* 
spirals have a constant V rot at large radii, whereas the curves of 
sub-L* galaxies are rising even beyond a characteristic radius, 
and very luminous objects have a negative gradient at large 
radii. The characteristic radii which define a characteristic rota- 
tion velocity that can be used for the TF analysis differ slightly 
between the two approaches. In the URC91 form, this radius is 
equal to 2.2 scale lengths, whereas it is as large as 3.2 rj for the 
URC96. Throughout the following sections, the quantity V max 
will refer to the usage of equation^ whereas V2.2 and V opt de- 
note an input of the URC91 and URC96, respectively. 

In the next step of the simulation, the two-dimensional ve- 
locity field was generated on the basis of the respective intrin- 
sic rotational law (e.g. Eq.^&EJ. tilted and rotated according 
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Fig. 4. Simulated velocity field for galaxy #4922 inclined by 
i = 64° with an intrinsic rotation curve of V max = 156 km/s. Thin 
solid and dotted lines indicate iso-velocity zones ranging from 
+ 120 km/s to -120 km/s rotation in line-of-sight, i.e. compris- 
ing a factor sin/. The thick solid lines denote the slit, thick 
dashed lines sketch the isophotes of the galaxy, with the outer- 
most corresponding to 2<x sky noise level in the 3000 s refer- 
ence /-band image (see Sect. 15. li . The misalignment between 
slit and apparent major axis is 6 = 13°. 



to the observed inclination and position angle of the respec- 
tive galaxy. The field was then weighted by the normalized sur- 
face brightness profile, i.e. brighter regions contribute stronger 
to the velocity shift at a given radius. Following this, the field 
was convolved. The Point Spread Function was assumed to be 
Gaussian with a FWHM determined from the mean DIMM val- 
ues during the spectroscopy. These values had to be slightly 
increased depending on the redshifted wavelength of the emis- 
sion line from which the observed RC was derived. We de- 
termined the correlations for this increase by measuring the 
FWHM on sets of VLT exposures in B, R, /, and comparing 
these to the according DIMM values. After the double folding, 
a strip of one arcsecond width was extracted from the velocity 
field at an angle that matched the slit misalignment in the obser- 
vation (see Fig.0J. The final computation step was an integra- 
tion perpendicular to the simulated slit, i.e. the projection of the 
velocity field strip onto the spatial axis. In all of the following, 
we will refer to this resulting position-velocity model as the 
synthetic rotation curve, whereas the rotational law (e.g. Eq.^ 
that is used as input for the modelling procedure will be re- 
ferred to by the term intrinsic rotation curve. The observed RC 
as described in Sect. !4.3l would directly reproduce the intrinsic 
RC only for a disk of 90° inclination (i.e., perfectly edge-on) 
that is observed with an infinitely thin slit at infinitely large 
spatial and spectral resolution. 
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In our approach, the intrinsic quantity V max is the only free 
parameter that tunes the reproduction of an observed RC by 
a synthetic curve. Only in the case of 3 spirals, the gas scale 
length had to be kept as a second free parameter, since the val- 
ues based on Eq.[2]were too large. For the complete data set, we 
derived V max by a visual comparison of synthetic and observed 
RC and, alternatively, via a ^-fitting procedure based on the 
errors from the RC extraction. These two methods are com- 
pared in Fig. |5] The sample of 77 spirals is subdivided accord- 
ing data quality: Curves which clearly probe the region of con- 
stant rotation velocity at large radii are considered high quality 
(M,t>j = 36), whereas RCs with smaller extent or asymmetries 
are included in the low quality sub-sample (N \,i = 41). 

The error on V max is assumed to be 



+ V l 



.(tani) 2 fff + V* 



(3) 



Here, the first term on the right hand side is the error from the 
;^ 2 -fits of the synthetic to the observed RCs, covering the range 
3 km/s < cr x i < 59km/s. The last two terms are the propagated 
errors of the respective uncertainties of the inclination and the 
misalignment angle. To derive the contributions of the errors 
cr, and cr (5 , we used a simple geometric correlation between the 
observed and intrinsic rotation velocity, 



V^ = V^ s (sinir 1 (cos<J)- 



(4) 



For the high quality sample, the absolute and relative errors on 
V max fall in the respective ranges 3 km/s < o> max < 135 km/s 
and 0.03 < crv m „/V max < 0.61, with a median 
<^v m »/V max ) = 0.19™ 

As can be deduced from Fig. |5J the by-eye "fits" and the 
^ 2 -fits are consistent within the errors for the majority of the 
objects. Nevertheless, a systematic trend towards low values 
of V maXvt ,2 with respect to the visually derived V max , v i s is evi- 
dent for slow rotators. An inspection of the observed RCs re- 
vealed that the discrepancies mainly arise in cases of asymmet- 
ric shapes. The^ 2 -fits also are weighted towards the inner parts 
of the curves by the higher S/N of the emission lines at smaller 
radii. This is a disadvantage since the outer parts of an observed 
RC are the most robust source of V max in the modelling proce- 
dure. Moreover, we found no correlation between the reduced 
X 2 values and our definition of RC quality, except for curves 
of perfect symmetry. For these reasons, we used the visually 
derived values for the TF analysis. A further discussion of this 
topic will follow in Sect. 18. 31 

In cases of multiple usable emission features in a spectrum, 
the RCs based on different lines were mostly consistent within 
the errors. For these objects, V max was derived from the curve 
with the largest covered radius and highest S/N. 

We show a consistency check of the three alternatives for 
the intrinsic RC shape in Figs. [6] and The results using the 
simple "rise-turnover-flat" shape via Eq. ^ are in agreement 
with the URC91 to within 5% for 79% of the FDF spirals and 
to within 10% for 95% of the galaxies, without a detectable 
dependence on the absolute values. The URC96 yields rotation 
velocities at R opt w 3.2 r& which on the mean are larger by 7% 
than V max , with slightly increasing differences towards slow ro- 
tation. This partly is an effect of the different characteristic radii 
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Fig. 5. Comparison between the V max derivation by visual 
alignment of observed and synthetic rotation curves (x axis) 
and via ^ 2 -fits (y axis). Typical error bars are shown for two 
objects. 



in the two parameterisations of the universal rotation curve, 
which correspond to 2.2 for the URC91 and 3.2 for the 
URC96, respectively. However, the slight differences between 
Vmax and V opt do not affect the results from our TF analysis, see 
Sect. 



5. Photometry 

5.1. Luminosity profiles 

Since the FDF imaging had been done under varying seeing 
conditions (ranging, e.g., from 0.46" to 0.89" FWHM in the 
/-band, see Heidt et al. 2003 for a description of the image 
stacking), only a limited number of images was used for the 
measurements of the structural parameters. To combine high 
spatial resolution with sufficient S/N at the outer isophotes of 
the TF objects, the 10 /-band frames with the best seeing were 
co-added. This yielded a reference image with 3000 s total ex- 
posure time, 0.49" FWHM and a 50% completeness limit of 
-25.1'". 

The disk light distributions were fitted with exponential 
profiles. A bulge component could not be accounted for at 
the given resolution. We considered slight variations of the 
PSF across the reference frame with six stars in the range 
18 < / < 20. Our algorithm minimizes x 1 in the parameter 
space span by inclination, position angle, scale length and cen- 
tral flux. The sample galaxies cover a range of 25° < i < 80° 
with a median of (i) = 53° and a median error of (cr,) = 5°. 

As a test of the accuracy of our ground-based luminosity 
profile analysis, we also performed measurements on a "VLT 
simulation" of the Hubble Deep Field North (HDF-N, Williams 
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Fig. 6. Comparison between rotation curve fitting with a sim- 
ple "rise-turnover-flat" shape (y axis) and the universal rota- 
tion curve shape as introduced by Persic & Salucci ( I1991I x 
axis). Solid symbols denote high quality curves which cover 
the region of constant rotation velocity at large radii. Typical 
error bars are shown for two objects. 



Fig. 7. Comparison between rotation curve fitting with a sim- 
ple "rise-turnover-flat" shape (y axis) and the universal rota- 
tion curve shape as introduced by Persic et al. (1996 x axis). 
High quality curves (solid symbols) cover the region of con- 
stant rotation velocity at large radii. Typical error bars are 
shown for two objects. 



et al. 119961 . The original drizzled images were re-binned to a 
scale of 0.2" per pixel and convolved with a Gaussian PSF of 
0.49" FWHM to match the characteristics of our reference I- 
band frame. We selected 40 objects with a variety of B/T ratios 
and i > 30° from the electronically available data published 
by Marleau & Simard dl998> . The authors did apply two- 
component fits to more than 500 galaxies in the HDF-N using 
the GIM2D package (Galaxy Image Two-Dimensional, Simard 
et al. l2002> . We found that for B/T < 0.4, the mean difference 
between the inclinations derived from the "VLT simulation" 
frame and the original GIM2D values was Ai = 1.8° + 7.7°. 

The median inclinations of our TF sample sub-divided ac- 
cording to SED type are <z) = 51° for T < 3, <i> = 53° for 
T = 5 and (i) = 51° for T — 8, i.e. potentially unresolved bulge 
components of early-type spirals in our ground-based imaging 
did not introduce a detectable bias. 

5.2. Rest-frame magnitudes 

Total apparent magnitudes were derived with the Source 
Extractor package (Bertin & Arnouts 1996 1. This program 
offers different algorithms for the photometry. We used the 
Magjxuto routine with variable elliptical apertures which is 
based on the "first moment" algorithm by Rron ( 1980 1 since 
it best reproduces the total magnitudes of extended sources. 

We will now briefly discuss the issue of the ^-correction. 
If available broad band information is limited to a few filters 
(in most previous studies two or even only one), the different 



wavelength ranges covered by a passband in observed frame 
and rest-frame introduce a strong dependence of the possible 
^-correction accuracy on SED type. At redshift z = 0.6, e.g., the 
difference between T — 4 and T — 6 corresponds to a change 
of half a magnitude in the transformation from B ODS to B rest , 
k B (e.g. Frei & Gunn 1994 1. Values of k B differ even more for 
earlier types, thus even a slight misclassification can introduce 
a substantial offset in the derived luminosity if observations are 
limited to one or two filters. In contrast to this, our TF project 
greatly benefits from the multi-band imaging of the FDF: The 
photometry in B, g, R and I enables us to use the filter that 
best matches the rest-frame fi-band to transform an apparent 
magnitude X into absolute magnitudes Mb up to the highest 
redshifts in the sample. 

We computed the ^-correction for our objects via synthetic 
photometry on SEDs in the range 1 < T < 10. As templates, we 
used the spectra published by Moller et al. (2001 1 which were 
generated with evolutionary synthesis models. The SEDs were 
redshifted by re-calculating the original flux Fq at wavelength 
A according to 



FoW(l+z)] 



(5) 



(1+z) 

(e.g. Contardo, Steinmetz & Fritze-v. Alvensleben 1998 1. 
Transformation from the apparent magnitude of a spectrum of 
type T at redshift z observed with a FORS filter X to the un- 
redshifted spectrum in Johnson B yields ks(X, T, z). Respective 
filters used for the input magnitudes were Bfors for z < 0.25, 
Sfors for 0.25 < z < 0.55, R F0RS for 0.55 < z < 0.85 and 
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^fors for z > 0.85. This ^-correction is much less sensitive 
to spectral type than a transformation B Q ^ — » B rest , especially 
at higher redshifts: e.g., k B (R, 1,0.6) - k B (R,5,0.6) = 0.05'" 
whereas k B (B, 1 , 0.6) - k B (B, 5, 0.6) = 0.64 m . 

For testing purposes, we additionally derived the colors 
of the templates purely within the Johnson-Cousins Filter 
system and compared them to the values published by 
Fukugita et al. d!995t . finding typical absolute deviations of 

0. 03 m < |A(X - Y)\ < 0.08" 1 . 

The second critical correction that has to be applied to the 
observed magnitudes is the inclination-dependent intrinsic ab- 
sorption A' B by the dust disks of the objects. Like, e.g., Vogt et 
al. (11996111997b and Milvang -Jensen et al. J2003> . we adopted 
the approach by Tully & Fouque J1985l >. It is based on geo- 
metric assumptions and usable for inclinations up to i = 80°, 

1. e. for our complete sample. The dust disk scale height is as- 
sumed to be half the scale height of the luminous disk, relative 
dust content is independent of mass or type, with an optical 
depth of t = 0.55. For the objects in our sample, A' B spans val- 
ues from 0.30 m at i = 25° to 0.96"' at i = 80°. We emphasize 
that the absorption for face-on disks is finite in this convention: 
A' B = 0.27™ for i — 0°. Our multi-band photometry ensures to 
compute the absorption for the same rest-frame wavelength in- 
terval at all the covered redshifts. 

For the galactic absorption at the coordinates of the FDF, 
we adopted the values which are given in Heidt et al. (U003 1, 
ranging from A| = 0.076"' to A] = 0.035"'. Let D A (z, H ) be 
the distance modulus (e.g. Peebles 1993 1 in the concordance 
cosmology, then the transformation from total apparent magni- 
tude mx to absolute magnitude M B is given by 

M B = m x - D A (z, H ) - k B (X, T, z) - A B - A% (6) 

If one assumes an SED classification error of AT = 2 (i.e., 
an Sc spiral could be misclassified as an Sb and vice versa), cor- 
responding to an uncertainty of the ^-correction of cr* < 0.1'" 
for all covered types and redshifts, the errors in absolute mag- 
nitude become 



where cr x is the random photometric error in the respective fil- 
ter and cr A ^ the uncertainty of the intrinsic absorption correc- 
tion via error propagation from the inclination error. Thanks 
to the very deep imaging, the random photometric errors are 
only 0.01"' on average and 0.03"' at maximum for the FDF spi- 
rals. Based on a comparison of our calibrations with archived 
FORS zero points, we estimate the systematic photometric er- 
rors to be < 0.0 l m ; these are neglected. The total errors <tm b of 
our complete sample fall into the range 0.07'" < ct Mb < 0.21"'. 

6. The data table 

We present the data of the spirals from our sample in Tabled 
This table is available in electronic form via anonymous ftp 
to cdsarc.u-strasbg.fr. The respective columns have the 
following meaning: 

ID — Original entry number from the FORS Deep Field 
photometric catalogue, see Heidt et al. (2003) for, e.g., the 



J2000 positions or aperture magnitudes in U, B, g, R, I, J and 
K. 

i — Disk inclination derived by minimizing x 2 of a two- 
dimensional exponential profile fit to the galaxy image in a 
coadded /-band frame with 0.49 arcsec FWHM (see Sect. 15. 1> . 
In some cases, our initial constraint of i > 40° had to be re- 
laxed. 

5 — Absolute misalignment angle between the MOS slit 
and the apparent major axis. Our initial constraint of 6 < 15° 
could not always be met during the construction of the setups. 

— Apparent disk scale length in arcseconds, derived via 
the same fits as i and 6. While tests confirmed that is the 
least critical input parameter in the V max derivation process, it 
is probably the one which is affected the strongest by the lim- 
itations of the ground-based imaging. Hence, we recommend 
not to use these values for applications like the Fundamental 
Plane of spiral galaxies. 

z — Spectroscopic redshift. 

T — SED type in the de Vaucouleurs scheme (see Sect. 14. 21 
for a description of the classification criteria). T — 1 accounts 
to Hubble type Sa, T = 3 to Sb, T = 5 to Sc, T = 8 to Sdm and 
T = 10 to Im. 

X — Total apparent magnitude derived with the Magjiuto 
algorithm of the Source Extractor package (Bertin & Arnouts 
1996 1 on the coadded FDF frames. Depending on the galaxies 
redshift, the filter X was chosen to best match the rest-frame 
B-band. For the least redshifted spirals with z < 0.25, this was 
X - B, for 0.25 < z < 0.55 we chose X = g, for 0.55 < z < 0.85 
we set X = R, and for the highest redshifts in our sample with 
z > 0.85 we used X = /. 

k B — K-correction for the transformation from the respec- 
tive filter X (see above) to rest-frame B as computed via syn- 
thetic photometry, see Sect. l5.2l for details. 

A g x — Galactic absorption in the respective filter X (see 
above) as given in Heidt et al. (2003 1. 

A' B — Intrinsic inclination-dependent dust absorption in 
rest-frame B following Tully & Fouque (1985 1 with the con- 
vention of a non-negligible extinction for face-on disks, i.e. 
A' B = 0.27'" at i = 0°. 

Da — Distance modulus in concordance cosmology with 
Q m = 0.3, Q A = 0.7 and H Q = TOkms" 1 Mpc" 1 . 

M B — Absolute B-band magnitudes computed via M B = 
X — k B — A s x -A' B - Da- The errors include the uncertainties in 
X, k B and A' B , see Sect. l5.2l for details. 

V m ax — Intrinsic maximum rotation velocity derived via 
synthetic velocity fields (see Sect. |4j assuming a linear rise of 
the rotation velocity at small galactocentric radii and a flat RC 
at large radii ("rise-turnover-flat" shape). The errors on V max 
were computed according to Eq.|3] 

B-R — Rest-frame color index, corrected for intrinsic ab- 
sorption. For galaxies with z < 0.25, apparent R magnitudes 
were transformed to rest-frame R, while apparent / magni- 
tudes were transformed to rest-frame R for all other objects, 
see Sect. l8.1l for details. 

note — Label indicating the RC quality. High quality 
curves ("H") extend well out to the region of constant rota- 
tion velocity at large radii. RCs of low quality ("L") have a 
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Table 1. The data of the 77 spirals from the FORS Deep Field Tully-Fisher sample. 
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0.52 


0.6535 


5 


21.78 


-0.58 


0.04 


0.41 


42.97 


-21.07±0.09 


169±36 


1.10 


H 


2533 


65 


34 


0.90 


0.3150 


5 


20.42 


+0.04 


0.06 
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0.04 


0.61 
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H 
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59 


12 
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3 
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0.08 


0.48 
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H 
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69 


2 
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1 
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+0.46 


0.06 
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41.65 
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249±53 
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L 
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25 


39 
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3 
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75 
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H 
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61 
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L 


4498 


40 


40 


0.33 


0.7827 


8 


22.60 


-0.38 


0.04 


0.34 


43.45 


-20.85±0.11 


95±58 


0.82 


L 


4657 


37 


5 


0.70 


0.2248 


5 


21.79 


+0.73 


0.08 


0.33 


40.24 


-19.58+0.08 


212+14 


0.91 


H 


4730 


32 


3 


0.90 


0.7820 


1 


20.94 


-0.23 


0.04 


0.31 


43.44 


-22.63±0.08 


430±84 


1.41 


L 


4806 


65 


4 


0.70 


0.2214 


5 


22.07 


+0.72 


0.08 


0.56 


40.21 


-19.50±0.08 


209±17 


0.72 


H 


4922 


64 


13 


0.71 


0.9731 


10 


21.77 


-1.06 


0.04 


0.54 


44.03 


-21.78+0.11 


156±27 


0.22 


H 


5022 


73 


10 


0.61 


0.3385 


10 


22.94 


-0.04 


0.06 


0.72 


41.26 


-19.06±0.14 


75±6 


0.51 


L 


5140 


50 


5 


0.60 


0.2738 


5 


22.94 


-0.08 


0.06 


0.40 


40.73 


-18.16+0.10 


125+39 


0.85 


L 


5286 


65 


1 


0.34 


0.3337 


8 


23.88 


+0.06 


0.06 


0.56 


41.22 


-18.02±0.18 


23±7 


0.64 


L 


5317 


53 


5 


0.59 


0.9745 


5 


21.39 


-0.92 


0.04 


0.42 


44.03 


-22.18+0.09 


236±36 


0.71 


L 



Table 1. — continued. 
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smaller radial extent and partly feature signatures of moderate 
kinematic perturbations like waves or asymmetries. 

7. The distant B-band Tully-Fisher relation 

To be able to derive the spectrophotometric and/or kinematic 
evolution of spirals at intermediate redshifts, it is crucial to 
carefully choose a data set of spirals at low z that can be used 
as reference. Consistency between distant and local data set in 
terms of the intrinsic absorption correction is one of the key is- 
sues. As stated in the introduction, a large number of local sam- 
ples have been constructed during the last decade, with kine- 
matic data based on radio observations or optical spectra. For 
our purposes, the term "local" refers to redshifts below z < 0.05 
corresponding to systematic velocities V sys < 15000km/s. In 
Ziegler et al. (2002), we selected a sample by Haynes et al. 
(1999 1 as a local reference. Here, we will use the data set 
of Pierce & Tully (FT9921 PT92 hereafter) instead. Our initial 
choice was basically motivated by the very good statistics of 
the Haynes et al. sample which comprises approx. 1200 spirals 
mainly of type Sc, which is also the most frequent SED type in 
our distant sample. 

On the other hand, using the PT92 data makes our results 
directly comparable to the studies of e.g., Vogt (|2001 1 and 
Milvang-Jensen et al. ©03 1. Moreover, the TFR of PT92 ben- 



efits from the inclusion of spirals with Cepheid-calibrated dis- 
tances. As for our sample, the photometry has been corrected 
for dust-reddening via the Tully & Fouque approach. A differ- 
ence that has to be accounted for is the convention for face- 
on extinction (see Sect. 15. 2> used by PT92, i.e., an offset of 
AMb = 0.27™ has to be applied. This way, the local fi-band 
TFR transforms into 



Mb — —7.48 log V max — 3.52 



(8) 



with an observed scatter of <x B = 0.41'". This is in good agree- 
ment with a bisector fit to the Haynes et al. sample restricted to 
1097 galaxies with H I profiles classified as good quality by the 
authors: 

Mb — —7.85 log V max — 2.78. (9) 
Here, we transformed apparent /-band magnitudes into B via 



colors from Frei & Gunn ( 1994 1. The two TFRs are perfectly 
consistent in the low-mass regime and show only a small offset 
of AM * 0.18'" at V max = 300 km/s. 

In Fig. [8] we show the TF diagram of the complete FDF 
spiral sample along with the local relation from PT92. As in 
Sect. |4] our sample is sub-divided according to the RC qual- 
ity. For the high quality data, the curves have a sufficient spa- 
tial extent to probe the region of constant rotation velocity at 
large radii, thereby yielding robust values of V max . In the case 



12 



A. Bohm et al.: The Tully-Fisher relation at intermediate redshift 




Fig. 8. FORS Deep Field sample of spirals in the range 0. 1 < 
Z < 1.0 in comparison to the local TFR by Pierce & Tully 
( 1992 dotted lines give 3 cr limits). The distant sample is sub- 
divided according to rotation curve quality: High quality curves 
(solid symbols) extend well out to the region of constant rota- 
tion velocity at large radii. 

of the low quality data, on the other hand, we cannot rule out 
the possibility that at least for a fraction of the objects V max is 
underestimated, i.e. offsets from the local TFR towards higher 
luminosities could be overestimated. For this reason, we will 
use only the high quality data for the following analysis. These 
spirals cover the ranges -18.6'" > M B > -22. 2 m in absolute 
magnitude and 74km/s < log V max < 295km/s in maximum 
rotation velocity. 

For a given V max , the difference in luminosity of an FDF 
object from the local PT92 fit is given by 

AM B = 7.48 log V max + 3.52 + M B . (10) 

These offsets are shown as a function of redshift in 
Fig. [5] Although the scatter of the offsets is reduced from 
cam = 1.32'" to cr AM = 0.97"' when restricting the FDF sample 
to high quality RCs, this is still over a factor of two larger than 
for the local data set. We speculate that this partly is an effect 
of the observational limitations for distant spirals like, e.g., the 
low intrinsic spatial resolution, but also reflects a broader range 
of star formation efficiencies than in the local universe. This 
interpretation is supported by the smaller scatter of the offsets 
which originates from the uncertainty in V max and amounts to 
0.63 mag for the HQ data. 

A linear^ 2 -fit to the high quality data yields 

AM B = -(1.22 + 0.56) -Z- (0.09 + 0.24). (11) 

For this fit, the error was computed as 

°AM S - °mF df + 7 - 482<T k> g yS F 



Fig. 9. Offsets of the FORS Deep Field sample from the local 
TFR as given by Pierce & Tully ( 1992 1 as a function of redshift. 
Filled symbols denote rotation curves which extend well out to 
the region of constant rotation velocity at large radii (labeled 
high quality data). 

+<r 2 <T +7.48V 2 ogy p Ii , (12) 

where the second and fourth term are the propagated errors 
from the uncertainties in log V max for the FDF spirals and the 
PT92 sample, respectively. 

As can be deduced from Eq. ^2 we observe an increas- 
ing brightening with rising look-back time. This is expected 
as an effect of the younger stellar populations, i.e. a higher 
fraction of high-luminosity stars than in the local universe. 
Our result is in agreement with those of Barden et al. (|5003 
AM B = -1.1 ± 0.5 at (z) * 0.9) and Milvang- Jensen j2"UU3l . 
who finds AM B = -(0.9 ± 0.3) • z, note that the latter correla- 
tion would be slightly steeper in the cosmology adopted here. 
A much smaller brightening of less than 0.2 mag at z < 1 is 
found by Vogt (I2U0TV 

Besides the dependency of AM B on redshift, the compar- 
ison of our sample with PT92 in Fig. [8] indicates a correla- 
tion between the TF offsets and maximum rotation velocity. 
Even restricting to rotation curves which probe the "fiat" re- 
gion at large radii, a number of distant spirals in the regime 
V max ~ 1 00 km/s are overluminous with >3 <x confidence, given 
the observed scatter of 0.41 mag for the local sample. This can 
be seen in Fig. EH where the offsets are plotted against the log- 
arithm of V raax . A linear ^ 2 -fit to the high quality subsample 
with an error estimation as defined in Ea.ll2lvields 

AM B = (4.40 + 0.69) log V max -(10.31 + 1.55), (13) 

corresponding to a brightening by more than two magnitudes 
for the least massive spirals in our sample and negligible offsets 
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Fig. 10. Offsets of the FORS Deep Field sample from the local 
TFR as given by Pierce & Tully ( 1992) as a function of max- 
imum rotation velocity. We observe large overluminosities for 
objects with slow rotation (i.e. low mass) and negligible offsets 
for fast rotators (i.e. high-mass spirals). Filled symbols denote 
rotation curves which extend well out to the region of constant 
rotation velocity at large radii (labeled high quality data). 



-i i | i i i r 

High quality RCs • 
Low quality RCs G 




0.2 0.4 0.6 0.8 1 

Fig. 11. Offsets of the FORS Deep Field sample from the lo- 
cal TFR as given by Pierce & Tully ( 1992), plotted against the 
relative errors in V max . The weak correlation between these two 
parameters can attributed to the lower RC quality and the larger 
relative errors of the RC extraction for slow rotators, see text 
for details. Filled symbols denote rotation curves which extend 
well out to the region of constant rotation velocity at large radii. 



at Vmax ~ 220 km/s, which on the basis of Eq. |F]is the typical 
rotation for local spirals of luminosity 2-3L*. 

8. Discussion 

Before interpreting the possible physical implications of 
Eq.^J we want to comment on a potential correlation between 
the TF offsets and the errors in V max . In Fig. ^2 the offsets are 
plotted against the relative errors <xy max /y max . Even if only the 
galaxies with high quality RCs are considered, the overlumi- 
nosities seem to be higher for objects with larger V max errors. 
We see two basic reasons for this slight dependency. 

Firstly, the RC quality on the mean is lower for low-mass 
objects. In particular, the least massive FDF spirals all are clas- 
sified as low quality data. A lower rotation curve quality in turn 
leads to higher values of ay> (Eq.[3) and an increased total error 
of the maximum rotation velocity. And secondly, the errors of 
the Gaussian fits to the emission lines which are performed in 
the process of the RC extraction do not depend on the magni- 
tude of the velocity shifts. On the other hand, for a given radius 
and fit uncertainty, the relative errors on V°^(r) are smaller for 
fast rotators. Since the observed relative errors contribute to the 
derived value of <t x i, this error is on the mean larger for slow 
rotators. 

The combination of these two effects leads to larger relative 
errors for small values of V max , i.e., for spirals with large TF 
offsets. 



We will now consider potential biases due to galaxy-galaxy 
interactions or sample incompleteness, the impact of the intrin- 
sic rotation curves shape and the question of different conven- 
tions for the intrinsic absorption correction. 

8.1. A bias due to environmental effects? 

To some extent, a correlation between rotation velocity and TF 
offsets is to be expected from previous studies. Kannappan et 
al. (2002) found a color-residuals relation that reflects over- 
luminosities of blue spirals and argued that this could be at- 
tributed to enhanced star formation. Since galaxies with blue 
colors, i.e. late types, feature on the mean lower values of V max 
than red, early-type spirals (see Sect. I8.5I . a correlation be- 
tween colors and TF offsets should coincide with a relation be- 
tween the offsets and V max . 

To look into this effect, we computed the B - R color in- 
dex of our FDF galaxies in rest-frame. In contrast to the initial 
derivation of the absolute magnitudes, we did not use total ap- 
parent magnitudes but brightnesses that were measured within 
apertures of two arcseconds diameter on coadded frames which 
were convolved to the same seeing (see Heidt et al. |2003 for 
details). Similarly to the procedure for the Z?-band luminosi- 
ties, we transformed the observed /?fors into ^cousins on ^ ^ or 
low-redshift galaxies (z < 0.25) and used the transformation 
^fors — > ^cousins instead f° r all objects at larger distances. The 
absorption coefficients with respect to rest-frame R were com- 
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Fig. 12. Offsets of the high quality FORS Deep Field sample 
from the local TFR by Verheijen (2001, circles), compared to 
the residuals of a free bisector fit to the high quality FDF data 
(triangles), both as a function of rest-frame color. The solid 
line represents the color-residual relation of the local Verheijen 
sample, the dotted lines denote the 3 <x limits. Typical error bars 
are shown for two FDF objects. 

puted following Cardelli et al. (1989 1 using again the intrinsic 
absorption convention by Tully & Fouque ( 1985 1. 

For the purpose of the color-residual relation, we compare 
our data to the local sample from Verheijen (2001). We show 
the offsets with respect to the Verheijen TFR 

Mb — —8.1 log V m ax — 2.07 (14) 

for our high quality sample in Fig. [21 Only a weak dependence 
of the offsets on rest-frame color is observed. In contrast to the 
results for the local spirals, offsets of the FDF spirals almost 
solely populate the regime of overluminosities with a median 
of -0.7"' and 12 out of 36 galaxies feature offsets of more than 
one magnitude from the local TFR. Note that this is consistent 
with the offsets from the PT92 relation which have a median of 
-0.77'". 

The situation is changed if one assumes that the TFR slope 
and zero point could vary with cosmic time, i.e. if a free fit is 
applied to the distant sample. On the basis of the high quality 
data, a bootstrap bisector fit with 100 iterations yields 



M B = -(4.66 + 0.67) log V max - (10.43 ± 1.50). 



(15) 



To be consistent with the Verheijen sample, we performed the 
same fit using the intrinsic absorption convention by Tully et 
al. (Ea.ll7> which yielded 



M B = -(5.22 + 0.69) log V max - (9.08 + 1 .56). 



(16) 



New residuals using this relation for the FDF sample were de- 
rived analogous to Eq. [TO] and are compared to the initial val- 
ues in Fig. [21 The residuals computed via Eq. [H[]are mainly 



distributed within the range -1 < AMb < +1 around zero, 
similar to the results given by Verheijen, though the scatter of 
the distant sample is larger and the correlation between color 
and residuals relatively weak. We conclude that the large off- 
sets we find using Eq.no| i.e. under the assumption of the local 
TFR slope and zero point, can hardly be attributed to the color- 
residual relation of local spirals. 



Alternatively to an evolution of the TFR with look-back 
time, one possible interpretation could be that a fraction of the 
distant spirals are subject to galaxy-galaxy interactions, which 
could result in TF offsets of up to several magnitudes as stated 
by Kannappan et al. (2003 1. To perform a search for galaxy pair 
candidates within our TF sample, we combined our complete 
data set (i.e. all galaxies with derived redshifts) with the lower 
resolution spectra from the FDF high-z campaign (Noll et al. 
2004 1, yielding a total of 267 galaxies at z < 1 . As upper limits 
on the three-dimensional separation of two pair candidates, we 
adopted the results of Lambas et al. (2003 1. Based on a data 
set of ~ 10 5 objects from the 2dF Galaxy Redshift Survey, the 
authors found that a projected distance D pro j < 100 kpc and a 
relative radial velocity A V sys < 250 km/s are reliable upper lim- 
its to select galaxy pairs with enhanced specific star formation. 



12 spirals from the TF sample show spectroscopically con- 
firmed neighbors within these limits. Of these galaxies, two 
are also possible members of the cluster at z — 0.33 lo- 
cated in the southwestern corner of the FDF. We show the 
TF offsets of these pair/cluster candidates in comparison to 
the rest of the sample in Fig. The two sub-samples seem 
to be similarly distributed. This is also implied by a two- 
dimensional Kolmogorov-Smirnov test which yields a propa- 
bility of 0.31 that both subsamples are drawn from the same 
distribution function. We only observe a slight overrepresenta- 
tion of pair/cluster candidates towards low rotation velocities. 
These are all spirals with RCs classified as low quality data, i.e. 
with relatively small radial extention and/or asymmetries. In to- 
tal, only 5 out of 18 (corresponding to 28%) of the pair/cluster 
candidates have high quality curves, whereas for the rest of the 
sample this fraction accounts to 53%. This is in agreement with 
the results of Kannappan et al. ( 2003 1, who find that a fraction 
of galaxies in close pairs show asymmetric or truncated RCs. 
We therefore reject the idea that a significant fraction of the 
galaxies with high quality RCs are subject to interactions. 



The nearest neighbor search cannot yield all pairs at z < 1 
within the FDF, since we do not have spectroscopic redshifts 
on all galaxies in this regime (both the TF and the high-z study 
have a limit in apparent brightness). However, the aim of this 
test was to clarify whether the candidates which are identified 
differ from the rest of the sample. Since the five pair/cluster 
candidates included in the analysis show only moderate TF off- 
sets, we conclude that it is unlikely that our high quality sub- 
sample introduces a bias towards overluminosities due to grav- 
itational interactions. 
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Fig. 13. Offsets of the FORS Deep Field spirals from the lo- 
cal TFR as given by Pierce & Tully ( 1993 as a function of 
maximum rotation velocity. The FDF sample is sub-divided 
according to environment: Squares denote potential members 
of a cluster at z = 0.33, triangles are objects featuring neighbor- 
ing galaxies within AV sys < 250 km/s and D plo j < 100 kpc and 
circles are probably isolated galaxies. Filled symbols denote 
objects with rotation curves that extend well into the region of 
constant rotation velocity at large radii. 



8.2. Potential incompleteness effects 

In the following, we will adress the question of whether the 
deficiency of our sample of objects with Mb < -19 and 
log V max > 2.1 (cf. Fig. [8} may be an observational bias. 

For a magnitude-limited TF sample, not all objects within 
the field-of-view that are geometrically suitable enter the fi- 
nal data set. Towards the faint end of the observed luminos- 
ity distribution, the sample therefore is incomplete. At a fixed 
maximum rotation velocity, this may give rise to a bias against 
low-luminosity spirals and, in turn, lead to an underestimated 
TFR slope as pointed out by Giovanelli et al. J1997t . 

Our apparent brightness limit of R < 23'" corresponds to a 
limit in luminosity that is higher for larger redshifts. An impact 
of the incompleteness bias on our results should therefore coin- 
cide with a decrease of the TFR slope with increasing redshift. 
However, if we split the complete FDF sample into sets for ob- 
jects with z < 0.45 (41 galaxies) and z > 0.45 (36 galaxies), the 
bootstrap fit slopes we find for the two are a = -3.75 + 0.44 
and a = -3.77 + 0.59, respectively. 

Potentially, our sample may introduce a bias against fast 
rotators of low luminosity. This is because the spatially re- 
solved emission lines of large spirals (which on the average 
have higher V max values) cover a larger CCD area and thus have 
lower signal-to-noise ratios per pixel than small spirals. For a 



given total line flux, the data set could therefore preferably con- 
tain slow rotators. 

Since late-type spirals on average have stronger emission 
lines (and lower V max values, see Sect. 18. 5> than early-type spi- 
rals, this effect would induce a type-dependency of the TFR 
slope. For SED types T < 3, T = 5 and T > 8, we find slopes 
of a = -3.74 + 1.40, a = -5.18 + 0.63 and -4.78 + 0.45 (the 
large uncertainty of the first of these fits shows the influence of 
low number statistics). Although the 14 spirals of type Sb or 
earlier have a flatter tilt than the two other sub-sets, we do not 
find a significant correlation between spectrophotometric type 
and TFR slope within the derived errors of the bootstrap fits. 

It is therefore unlikely that the shallower slope of the FDF 
sample with respect to the local data is an observational effect. 
We emphasize that, if the observed lack of distant spirals in the 
regime of low luminosities and fast rotation was due to such a 
bias, then the intrinsic scatter of the distant TF relation would 
have to be increased by a significant factor compared to the 
local universe, which we do not see. 

8.3. Impact of the intrinsic RC shape 

As stated in Sect. 14.41 the rotation velocities derived via the 
RC modelling are consistent for the "rise-turnover-flat" shape 
and the universal rotation curve from Persic & Salucci ( 1991 1, 
but slightly different if the URC as given in Persic et al. {1996 1 
is used. However, if we restrict our sample to the high qual- 
ity data, the difference between V max and V op t is negligible 
for spirals with V max > 150km/s, and amounts to only 5% 
in the median for the slow rotators. In effect, the TF offsets 
would be altered towards lower luminosities by only ~0.15 m 
at V pt = 100 km/s if the URC96 was used alternatively. None 
of our results would be affected significantly by such a small 
difference. For this reason, our analysis will be based on V max 
throughout the section. 

Another topic that has been referred to earlier are the slight 
differences between ^ 2 -fits and "by-eye" fits of the synthetic to 
the observed rotation curves. Since the former yielded system- 
atically lower maximum rotation velocities than the latter, the 
TF offsets would be brighter by ~ 0.25'" at V max = 100 km/s 
if the analysis was based on the V max values from ;f 2 -fit RC 
modelling as introduced in Sect. 14.41 Therefore, if one assumes 
a "rise-turnover-flat" shape of the intrinsic RC, our approach 
of a by-eye comparsion between synthetic and observed RC 
yields conservative values of the TF offsets. 

8.4. Influence of the intrinsic absorption correction 

The convention we use for intrinsic absorption is purely 
inclination-dependent, as the optical depth and fractional dust 
disk thickness are held fixed for the entire sample. More re- 
cently, Tully et al. ( 1998 1 and Karachentsev et al. (2002) have 
found evidence for an internal extinction law which also de- 
pends on rotation velocity. These observations favour a higher 
amount of absorption in fast rotators than for spirals of low 
Vmax- If one assumes the relation 2V max ~ W' R between max- 
imum rotation velocity and H I profile linewidth (see Tully & 
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Fouque 1985 1, then equation 1 1 given in Tully et al. transforms 
into 



-6 



A V,i 
A B 



2.75(logV max -1.63)log( fl /£), 



(17) 



where a and b are the apparent major and minor axes, respec- 
tively. E.g., for a highly inclined disk with i = 80°, this yields 
0.59'" at V max = lOOkm/s and 1.35"' at V max = 300km/s, 
whereas the initial Tully & Fouque approach gives a value of 
A l B = 0.96 m independent of the maximum rotation velocity. 

In effect, the slope of any TFR would be steeper if the in- 
trinsic absorption is accounted for based on Eq. El But, since 
this mass-depedency is linear in log V max /MB-space, the slope 
change would affect distant and local sample consistently. This 
can be quantitatively verified with Eqs.^](the local Verheijen 
TFR) and^](the FDF high quality sample) which are both de- 
rived using the Tully et al. convention. Neither the median of 
the FDF spiral offsets from the local TFR (0.7 mag in the Tully 
et al. convention vs. 0.77 mag in the Tully & Fouque conven- 
tion) nor the evidence for a change in the distant TFR slope 
(>3cr confidence level in both conventions) do significantly dif- 
fer between the two approaches. 

Since our results are therefore independent of the intrinsic 
absorption convention, we used the Tully & Fouque approach 
throughout the paper for the sake of direct comparability with 
the previous studies. 

8.5. A mass-dependent luminosity evolution? 

It is well known that in the local universe, a dependency of 
V max on galaxy type is observed. Blue, late-type spirals are on 
average slower rotators than red, early-type spirals (e.g. Rubin 
et al. I1985> . For our sample, sub-divided according to SED 
type, we find respective median values of (V max ) ~ 248km/s 
for T < 3, <V max ) ~ 145 km/s for T = 5 and <V max > * 89km/s 
for T > 8. Since we observe a correlation between the TF off- 
sets and V max , these V max distributions imply a dependency of 
the mean luminosity offsets on SED type as is illustrated in 

Fig.m 

A consequence of this is a potential selection effect for 
small samples which mainly comprise a certain sub-type. If 
a data set was biased towards late-type spirals due to the target 
selection on blue colors (like, e.g., Simard & Pritchet lT998> or 
strong emission lines (e.g. Rix et al. 1997), a considerable evo- 
lution in luminosity would be derived. On the other hand, if a 
data set preferably contains early-type spirals, i.e. large disks 
(Vogt et al. 1996 1, only a modest luminosity offset from the 
local TFR will be observed. 

A straightforward interpretation of the correlation between 
the luminosity offsets and the maximum rotation velocity we 
find could be a change of the TFR slope with look-back time. 
As given in Sect. 18.11 a bootstrap bisector fit to the high qual- 
ity FDF data yields a slope of -4.66 + 0.67, corresponding to 
>99% confidence for a shallower slope at intermediate redshift 
with respect to the local PT92 sample. To verify this with a 
different fitting method, we estimate the errors of the absolute 
magnitudes to be 



-3 







HQ d 


1 1 

ata, Sa Sb 


■ 






HQ data. Sc 








HQ d 


ata, Sd...Im 


▲ 




\. O 




Lo HQ data 




- 
- 


& "i 






- 
- 




• 








- 




o 


\* 










m 

O 

• r 


□ " 




1,1,1 




□ 


□ 



1.5 



2.5 



log V max [km/s] 



Fig. 14. Offsets of the FORS Deep Field sample from the local 
TFR as given by Pierce & Tully ( 1992) as a function of maxi- 
mum rotation velocity, sub-divided according to SED type into 
Sb or earlier (squares), Sc (cicrles) and Sd or later (triangles), 
respectively. Filled symbols denote rotation curves which ex- 
tend well out to the region of constant rotation velocity at large 
radii (labeled high quality data). 

The first term on the right hand side corresponds to Eq.^and 
the second term is the error propagation from the errors in V max 
based on the PT92 slope. A linear x 2 -fit to the high quality 
FDF then gives 



M B = -(3.55 + 0.41) log V max - (12.84 + 0.92), 



(19) 



i.e. the significance for a TFR slope change is even higher than 
on the basis of a bootstrap bisector fit. 

As an explanation for the shallower slope at redshift 
Z ~ 0.5, we will now discuss two alternatives. Firstly, our result 
may point to a luminosity evolution that depends on the max- 
imum rotation velocity. According to simulations by van den 
Bosch ( 2002 ), the total mass of a spiral galaxy within the virial 
radius can be estimated via 



"■"= 254 10 'S^)(^) 



(20) 



i.e. the correlation AMg oc log V max implies AM# oc logM v i r . 
For the spirals with high quality RCs, covering the range 
3.2 • 1O 1O M < M vir < 1.6 • 1O I2 M «M vir ) = 2.3 • 1O U M ), a 
linear ^ 2 -fit yields 



AM B = (1.14 + 0.27) log — -(13.46 + 3.09). 

M 



(21) 



°"tfr 



4+7.48V 2 ogW 



(18) 



Although the ground-based luminosity profile fits may place 
only upper limits on the scale lengths for some of the appar- 
ently smallest galaxies, we find evidence for a mass-dependent 
luminosity evolution which accounts to ~ 2"' in rest-frame B 
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for the least massive objects and is negligible for high-mass 
spirals. This implies that the redshift dependency we observe 
(Eq. 1111 is most probably a lower limit, since low-mass spi- 
rals — which show strong evolution in luminosity according 
to Eq. — at higher redshifts will fall beyond our apparent 
brightness limit. For the same reason, we cannot speculate on 
a potential evolution (more precisely, a decrease) of the mean 
galaxy masses with redshift. 

If one assumes that low-mass spirals have not undergone 
a significant increase of their metallicities since z ~ 0.5, 
our result is similar to that of Kobulnicky et al. (|2003l, 
who used the distant luminosity-metallicity relation and found 
AM B as -1. . . -2 mag. On the other hand, Vogt J200U does not 
find a TFR slope change with a TF sample from the same sur- 
vey (DEEP Groth Strip Survey). 

Besides a mass-dependent luminosity evolution, a second 
possible explanation for the flatter tilt of the distant TFR could 
be a strongly starforming galaxy population within our sample 
that contributes less to the local luminosity density, i.e. either 
these galaxies could be overnumerous or overluminous at in- 
termediate redshift. In terms of the luminosity function, this 
topic is well known as the "faint blue galaxy excess" (see Ellis 
1 19971 for an overview). Based on spectroscopic data from the 
Canada-France-Redshift-Survey, Driver et al. ( 1996 1 found that 
dwarfs may have faded by over one magnitude between the 
regime 0.2 < z < 0.5 and the local universe. These observa- 
tional findings could be understood in terms of a reionization 
era between redshifts z ~ 6 and z. ~ 1 which suppresses star 
formation in low-mass dark halos (M vir < 1O 9 M , Babul & 
Rees fT992l . 

However, only a small fraction of the high qual- 
ity FDF sample could be members of such a blue 
dwarf population since the former covers the range 
74km/s < log V max < 295km/s and the lower virial 
mass limit of 3 ■ 1O 1O M also exceeds the dark halo mass 
range of dwarfs. Furthermore, the star formation rates 
derived from [Oil] equivalent widths fall in the range 
O.SMoyr 1 < SFR < 18.3 M yr _1 , with a median of 
~4.5M Q yr~', typical for massive local spirals (e.g. Kennicutt 

If our result is interpreted in terms of a decreasing luminos- 
ity of low-mass spirals over the past ~ 5 Gyrs (which accounts 
to the look-back time at z ~ 0.5), it could be due to a mass- 
dependent evolution of the mass-to-light ratio, or even a mass- 
age relation. A direct comparison to predictions of stellar pop- 
ulation models is difficult since any evolution of the TFR intro- 
duces several competing effects. On the one hand, younger stel- 
lar populations coincide with a decrease of the mass-to-light 
ratio. On the other hand, since less gas has been consumed via 
star formation, the gas mass fraction most probably increases 
with look-back time, thereby increasing the mass-to-light ra- 
tio. Moreover, within the framework of hierarchical merging, 
disk sizes should be smaller for a given maximum rotation ve- 
locity and masses on the mean lower towards higher redshift 
(e.g. Mao et al. 1998 1. Based on the rotation velocity-size rela- 
tion, i.e. the correlation between V max and ra, our sample shows 
some evidence for slightly smaller disks at z ~ 0.5 (see Bohm et 



al. 2003). However, due to the limitations of the ground-based 
imaging, the significance of this evolution is low. 

In effect, a decrease of the disk sizes and an increase of the 
gas mass fraction would tend to shift distant spirals to the low- 
luminosity side of the local TFR, whereas a lower stellar mass- 
to-light ratio would result in a shift to the high-luminosity side. 
A domination of the first two processes for fast rotators could 
explain why a fraction of the high-mass FDF spirals are under- 
luminous in the TF diagram. 

As already stated, the correlation between luminosity evo- 
lution and redshift given in Eq.^Jis probably only a conserva- 
tive estimate. Because the disk sizes and gas fractions may also 
evolve, the pure luminosity evolution possibly is larger than 
approx. one magnitude between z = 1 and the local universe. 
E.g., chemically consistent evolutionary synthesis models by 
Moller et al. (2001 ) yield a luminosity evolution of ~ 1.5 mag 
in rest-frame B for an Sc spiral over this redshift range. 

Although the mild evolution of the scale lengths that might 
be deduced from the FDF sample is in compliance with the 
"bottom-up" scheme of structure growth within the Cold Dark 
Matter hierarchical model, a luminosity evolution that is larger 
for low-mass galaxies would be in contradiction to it. If our 
result is interpreted in terms of younger ages (lower forma- 
tion redshifts) for lower masses, the contradiction will be even 
more obvious. This is simply because small Dark Matter Halos 
should have formed earlier than large ones, and thus the stel- 
lar populations in galaxies of low mass should be older than 
those of high-mass systems. Semi-analytic simulations which 
account for a mass-dependent evolution of the stellar popula- 
tions also yield an increase of TFR slope with look-back time 
(e.g. Ferreras & Silk l200TV 

On the other hand, it is well established and also valid for 
our sample that the colors of spirals tend to become redder with 
mass (beginning of this section) which is not reproduced by 
simulations within the hierarchical scenario unless the spec- 
trophotometric properties of local spirals are used as a calibra- 
tion (e.g. Bell et al. l2003l . 

We therefore conclude that, if the TFR slope of our distant 
sample is related to a mass-dependant luminosity evolution, 
this would be at variance with the hierarchical merging sce- 
nario on small scales. 

9. Conclusions 

Using imaging data and spectroscopy taken with the ESO Very 
Large Telescope, we have derived structural parameters and re- 
solved rotation curves of a magnitude-limited sample of 77 
spiral galaxies in the FORS Deep Field. The objects cover the 
redshift range 0.1 < z < 1.0 and comprise all types from Sa to 
very late-type. Via a rotation curve modelling that takes into 
account geometric effects as well as seeing and optical beam 
smearing, the maximum rotation velocities have been derived 
and the distant B-band Tully-Fisher relation was constructed. 

We find evidence for a luminosity evolution with look-back 
time which amounts to a brightening of AMb ~ -1"' at redshift 
unity. Moreover, we observe a correlation between the lumi- 
nosity evolution and the total masses. The distant low-mass 
spirals are brighter by up to two magnitudes than their local 
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counterparts, whereas the luminosity evolution of high-mass 
systems is negligible. In effect, the slope of the Tully-Fisher 
relation at intermediate redshift is shallower than for local sam- 
ples. This may partly be caused by a population of small, star- 
forming galaxies that contribute less to the luminosity density 
in the local universe. Nevertheless, the vast majority of our ob- 
jects have virial masses too large to be dwarf galaxies and show 
star formation rates typical of normal spirals. The flatter tilt we 
find for the distant Tully-Fisher relation is in contradiction to 
the predictions of recent semi-analytic simulations. 
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